Z4— U(l) crossover of the order parameter symmetry 
in a two-dimensional valence-bond-solid 
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We discuss ground-state projector simulations of a modified two-dimensional S = 1/2 Heisenberg 
model in the valence bonds basis. Tuning matrix elements corresponding to the diagonal and off- 
diagonal terms in the quantum dimer model, we show that there is a quantum phase transition from 
the antiferromagnet into a columnar valence- bond- solid (VBS). There are no signs of discontinu- 
ities, suggesting a continuous or very weakly first-order transition. The Z4-symmetric VBS order 
parameter exhibits an emergent U(f) symmetry as the phase transition is approached. We extract 
the associated length-scale governing the U(l)-Z4 cross-over inside the VBS phase. 
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A valence-bond-solid (VBS) is a magnetically disor- 
dered state of a quantum spin system in which transla- 
tional symmetry is spontaneously broken due to the for- 
mation of a pattern of strong and weak bond correlations 
(Si • Sj) (where i,j are nearest-neighbor sites). Using 
an SU(A^) generalization of the Heisenberg model, Read 
and Sachdev showed that a four-fold degenerate colum- 
nar VBS ground state can be expected on the square 
lattice [f]. Numerical studies have found evidence for 
such VBS states in frustrated SU(2) symmetric systems 
, but because of technical limitations, in particular the 
sign problem in quantum Monte Carlo (QMC) simula- 
tions 0, the nature of the strongly- frustrated ground 
state remains controversial 0| • Another challenging issue 
is how the ground state evolves from an antiferromagnet 
(AF) into a VBS. According to the "Landau rules", one 
would expect a direct transition between these states to 
be first-order Q , because unrelated symmetries are bro- 
ken. There could also be an intervening disordered (spin 
liquid) phase Q or a coexistence region. Senthil et al. re- 
cently suggested an alternative scenario for a generic con- 
tinuous transition based on a "deconfined" quantum crit- 
icalpoint (DQCP) associated with spinon deconfincment 
0, Q- This proposal has generated significant interest, 
as well as controversy. An extended "J-Q" Heisenberg 
model has been introduced Q which is not frustrated, in 
the standard sense, but includes a four-spin interaction 
which destroys the AF order and leads to a VBS ground 
state. This model is amenable to large scale QMC stud- 
ies, which show scaling behavior consistent with a DQCP 

Other studies dispute these findings, however 
Numerical studies of the proposed field theory describing 
the deconfined quantum critical point are also subject 
to conflicting interpretations 1^, ll3| . Further studies of 
AF-VBS transitions is thus called for. 

In this Letter we address an important aspect of the 
VBS state and the AF-VBS transition, namely, the na- 
ture of the quantum fluctuations of the VBS order pa- 
rameter. In the DQCP theory, the Z4 symmetric lattice- 
imposed structure of the VBS is a dangerously irrele- 



vant, and, as a consequence, U(I) symmetry emerges 
close to the DQCP 0. An U(I) symmetric VBS order 
parameter was indeed confirmed in the studies of the J-Q 
model [^, [13, Hi, and also in simulations of the SU(iV) 
Heisenberg model with iV > 4 14]. However, the ex- 
pected cross-over into a Z4 symmetric distribution inside 
the VBS phase was not observed. This can be inter- 
preted as the lattice sizes studied so far being smaller 
than the spinon confinement length-scale A, which gov- 
erns the U(l)-Z4 cross-over 15]. A should diverge as 
where is the dimer (VBS) correlation length and a > I 
'T], and for a finite lattice with L ^ A the distribution 
should be U(l) symmetric. The models studied so far 
have a rather weak VBS order, and hence is large, 
which likely makes it difficult to satisfy L ^ The 
exponent a is not known. 

Here we introduce a way to generate much more robust 
VBS states, with which we can study the U(l)-Z4 cross- 
over already on small lattices. Our approach is based 
on a ground-state projector QMC method operating in 
the valence bond (VB) basis [3, 12, l3]- Starting from 
some trial state I^E*), the ground state of a hamiltonian H 
can be obtained by applying a high power of H; Y^ii) ^ 
H"^Y^). Consider the S — ^ Heisenberg model written 
as a sum of singlet projection operators -ffy. 
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where («, j) denotes nearest neighbors on a square lattice 
oiN = L? sites. In the VB basis the trial state \ '^) is a su- 
perposition of singlet products |(ai,6i) • • • (aAr/2, &Ar/2)), 
where {a,b) = {]aib — ia]b)/V^ with a and b sites on 
different sublattices. We here use the amplitude-product 
state of Liang et al. 1^, 2^. A singlet projector can have 
two different effects upon acting on a VB state; 
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These rules form the basis of the VB projector method 
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FIG. 1: (Color on-line) Diagonal (a) and oflt-diagonal (b) sin- 
glet projection operations (indicated by the arches) on VB 
pairs on a plaquette. The matrix elements ((2)| and ((3)l cor- 
responding to these operations are multiplied by Q„ and Qk, 
respectively. In (a) , the factor is 2Q„ if there is a VB also on 
to the left side of the operator. For all other bond configura- 
tions the matrix elements remain those in and Q. 



171 Il8l | , where i?™ is expanded in its strings of m sin- 
glet projectors. Expectation values of operators O are 
obtained by importance-sampling the VBs and operator 
strings produced when expanding 
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For details of these procedures we refer to Refs 

In Ref. 0], the J-Q model, which includes a four-spin 
coupling consisting of terms —QHijHki, with ij and kl 
site pairs on two opposite edges of a plauqtte, was stud- 
ied using the VB projector method. The Q term nat- 
urally favors singlet formation on plaquettes, and this 
was shown to lead to a VBS state when J/Q < {J/Q)c, 
{J/Q)c ~ 0.04. Here we introduce another mechanism 
leading to VBS formation. We define an effective hamil- 
tonian based on the Heisenberg model in the VB basis, by 
changing the diagonal matrix element 1 in Eq. ([2]) and the 
off-diagonal matrix element i in Eq. ([3]) to and ^Qfc, 
respectively, for Hab acting on VBs on opposite edges 
of the same plaquette. These operations, illustrated in 
Fig. [U correspond to the kinetic- and potential-energy 
terms of the quantum dimer model [21|. There, how- 
ever, the Hilbert space consists of only dimers connect- 
ing nearest-neighbor sites, whereas we here keep the full 
space of VBs connecting any pair of sites on different sub- 
lattices. In the quantum dimer model, the dimer configu- 
rations are also considered as orthogonal states, whereas 
we here keep the singlet nature of the VBs, whence the 
states are non-orthogonal. The non-orthogonality may at 
first sight seem problematic, because when Qv,Qk 7^ 1 
the hamiltonian is non-hermitean. We therefore refer to 
it as an pseudo hamiltonian in the VB basis. However, 
in spite of this, the states generated by the projection 
procedure (with the sampling weights modified by the 
presence of the factors Q„ and Qk, and taking the power 
m large enough for convergence to the m — oo limit) are 
completely well-defined SU(2) invariant quantum states. 
We can thus think of the modified projection technique as 
a means of generating a family of states parametrized by 
Qy and Qk- Moreover, there must be some corresponding 




FIG. 2: (Color on-line) Finite-size scaling of the spin and 
dimer correlation lengths. The inset shows ^/L versus the 
coupling, with crossing points tending toward Qy ~ 1.40. 

hamiltonians, defined in terms of the standard spin oper- 
ators Si, which have these states as their ground states. 
Although we are not able to write down these hamiltoni- 
ans (which likely contain multi-spin interactions, possibly 
long-ranged), it is still useful to study the evolution of the 
states as a function of Qv and Qk- Here we will consider 
two cases; Qv > l,Qk = ^ and Qk > l,Qv = 1, which 
we refer to as the Qv and Qk models, respectively. Both 
these models indeed undergo AF-VBS transitions. 

We calculate the square of the staggered magnetiza- 
tion, = (M • M), where 

M-^E(-1)'^'S.>^ (5) 

is the operator for the AF order parameter. The colum- 
nar VBS operator for x-oriented bonds is 



Dx — ^ ^^( ^)^Sx,y ■ Sx+l,yj 

x,y 



(6) 



and Dy is defined analogously. We calculate the squared 
order parameter, = {D^ + -^y)' distribution 
P{Dx,Dy) as in 0. Results for these quantities and 
the corresponding spin and dimer correlation lengths 
(spin) and (defined through the momentum-space sec- 
ond moments of the spin and dimer correlation functions) 
indicate coinciding critical points for the AF and VBS or- 
der parameters. In the following we first discuss the finite 
size scaling behavior of the Qv model. 

We define a reduced coupling q = Qv — Qv- Then, if 
there indeed is a single critical point, there is AF order 
for q <{) and VBS order for g > 0, and in the thermody- 
namic limit the squared spin and dimer order parameters 
should scale as ~ {—q)'^^^ and D"^ ^ q^^"^ inside the 
respective phases. To extract and the exponents, we 
use standard finite-size scaling forms, 



L'''il + aL-^)Fs{qL^^n, 
L'''{l + aL-^)Fd{qL'^n, 
L{l + aL~^)G,AqL^^l- 



(7) 
(8) 
(9) 
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FIG. 3: (Color on-line) Finite size scaling of the squared AF 
(top panel) and VBS (bottom panel) order parameters of the 
Q„ model. The insets show the unsealed data. 

where — 213s/ f^d — "^Pa/i^ and the correlation length 
exponent v is the same for all the quantities (as required 
in the DQCP theory). The scaling functions Fg^d and 
Gs.d are extracted in the standard way by adjusting the 
critical point and exponents to collapse finite-size data 
onto common curves. Since our lattices are not very 
large, L < 24, a subleading correction helps significantly 
to scale the data. In all cases we find that uj = 1 works 
well (the prefactor a is quantity-dependent, however). 

Results are shown in Fig. [2] and [31 All the data can 
be scaled with Ql = 1.400(5), v = 0.78(3), = 0.27(2) 
and /3d = 0.68(3). Here v and jid are approximately the 
same, within error bars, as in the J-Q model @, while /3s 
is very different — for the J-Q model /9s ~ /Sd = 0.63(2) 
was found. The range of system sizes is quite small and 
we cannot, of course, exclude drifts in the exponents for 
larger lattices, nor a very weakly first-order transition. 

Turning to the Qfe model, it is more demanding com- 
putationally, because the critical point is rather large, 
Qfe = 2.5(1), leading to a lower acceptance rate in simu- 
lations close to the transition than for the Q„ model. We 
can therefore not reach the same level of precision for the 
exponents. The results are nevertheless consistent with 
a continuous transition and exponents similar to those of 
the model. 

Unfortunately, we cannot easily calculate the dynamic 
exponent z with the present approach, because it requires 
access to the triplet sector, e.g., to extract the spin gap 
A ~ L~^. Our model is explicitly defined only in the 
singlet sector. While one can extend the VB basis and 
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FIG. 4: (Color on-line) VBS order=parameter distributions. 
The left column is for Qfe = 1, Q„ = 1.44, 1.48, 1.54 (from top) 
on 12 X 12 lattices. The right column is for = l,Qk = 
2.5, 3.5, 5.0 (from top) on 16 x 16 lattices. 



the projection scheme to triplets 1^, 17 1, the extension 
of the Qv and Qk models to this sector is not unique, and 
z may depend on how that is accomplished. We could in 
principle calculate gaps in the singlet sector, but this is 
much more complicated. 

Our main interest in studying the Qy and Qk models is 
in the distribution P(Dx, Dy) of the columnar dimer or- 
der parameter. While this is a basis dependent quantity, 
it still provides direct information on the order parameter 
symmetry. In a columnar symmetry-broken VBS state, 
we expect a distribution with a single peak located on the 
X OT y axis, while in a plaquette state the peak should be 
on one of the 45° rotated axes. In simulations that do not 
break the symmetry, we expect four-fold symmetric dis- 
tributions, with peak locations corresponding to the type 
of VBS as above. In previous studies of VBS states, only 
ring-shaped distributions were observed 0, [I^, [ll|, [l^ . 
however, which can be taken as a confirmation of the 
predicted emergent U{1) symmetry close to a DQCP. 
One would then expect the four-fold symmetry to appear 
for large systems, L ^ A, inside the VBS phase, as has 
been observed explicitly in a classical XY model includ- 
ing dangerously irrelevant Zq {q > 4) perturbations [3]. 
With the Qy and Qk models, we can reach further inside 
the VBS phases than in the previously studied quantum 
spin systems, and, as seen in Fig. 31 we can indeed follow 
the evolution from U{1) to symmetric distributions 
as a function of the coupling constants even for modest 
system sizes. The peak locations correspond to columnar 
VBS states for both models, although the shapes of the 
distributions are qualitative different in other respects. 
The previously observed ring-shaped distributions [9|, [14 1 
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FIG. 5: (Color on-line) Finite-size scaling of anisotropy order 
parameter. The inset shows the unsealed data. 



and the VBS order-parameter exponent j3d are roughly 
consistent with those found previously for the J-Q model 
(which is the best candidate so far for a DQCP), 
the AF exponent Ps is significantly smaller. It thus ap- 
pears that these transitions are in different universal- 
ity classes. This implies that emergent U(l) symmetry 
(which is associated with a dangerously irrelevant lattice- 
imposed potential) is a salient feature of VBS states that 
is more generic than the particular DQCP scenario by 
Senthil al. Our results then also point to a broader 
range of continuous (or, possibly, very weakly first-order) 
AF-VBS transitions. 

We would like to thank Naoki Kawashima and Masaki 
Oshikawa for stimulating discussions. This work was sup- 
ported by the NSF under grant No. DMR-0803510. 



are more reminicant of those for the Qfe model. 

To study the length scale A which governs the Z4- 
U{1) crossover (and is related to the scaling dimension 
of a dangerously irrelevant field [3|) we define an order [2 
parameter which is sensitive to the Z4 anisotropy 



D4 = 



J dD, J dDyP{D,,Dy)r,y cos{W) [3 



1 I-2TV 

dr / der^P{r,9)cos{Ae), (10) 



where r^y — {D^ + Dy)^/^. This order parameter should 

(11) 



! order p; 
obey the finite-size scaling form 

D4 ^ L''-''''^Fi{qL^/''^''), 



with 04 > 1. The data can be scaled with a = 1.30(5), 
as shown in Fig. [5] (where we use the same Q^^^ad, and 
1/ as in Figs. [3] and [2]). Here the error bars on the raw 
data, as seen in the inset, are much larger than for Z?^, 
reflecting that slow angular fluctuations of the VBS order 
parameter in the simulations (which do not affect the 
rotationally-invariant D^). In the classical XY model 
with Z4 perturbation w 1.1 [isl . [23j . and, thus, the 
Qy model has more prominent angular fluctuations. For 
the Qk model, we find an even larger 04 ~ 1.5, although 
the error bars are rather large and we cannot say for sure 
that it is different from the Qy model. 

To summarize, by tuning specific matrix elements in 
valence-bond QMC simulations, we are able to study a 
family of SU(2) symmetric states undergoing AF-VBS 
phase transitions. Unlike previous studies of quantum 
spin models with VBS states, we are able to observe 
both the Z4 symmetry of the order-parameter distribu- 
tion (which arises from the nature of the VBS on the 
square lattice) deep inside the VBS phase and the cross- 
over into an emergent U{1) symmetry upon approach- 
ing the transition point. We extracted the length scale 
A, which is associated with spinon deconflnement in the 
DQCP theory. While the correlation length exponent 
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